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We explore simple but novel bouncing solutions of general relativity that avoid singularities. These 
solutions require curvature k — +1, and are supported by a negative cosmological term and matter 
with — 1 < w < —1/3. In the case of moderate bounces (where the ratio of the maximal scale factor 
a+ to the minimal scale factor a_ is 0(1)), the solutions are shown to be classically stable and cycle 
through an infinite set of bounces. For more extreme cases with large a+/a_, the solutions can 
still oscillate many times before classical instabilities take them out of the regime of validity of our 
approximations. In this regime, quantum particle production also leads eventually to a departure 
from the realm of validity of semiclassical general relativity, likely yielding a singular crunch. We 
briefly discuss possible applications of these models to realistic cosmology. 
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Two questions have recurred again and again in theo- 
retical cosmology, starting with [1 -3J: 1) is the Universe 
eternal, or did it have a beginning at some definite time 
in the past?, and 2) is it possible to make Universes which 
enjoy one or more "bounces" where the scale factor first 
crunches, and then bangs? [T8] . 

The answers to these two questions are deeply inter- 
twined with the subject matter of the singularity theo- 
rems of Penrose and Hawking (discussed comprehensively 
in [4]). These theorems show that, given an energy con- 
dition of the form 
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for a suitable class of vectors w M , where T M „ is the stress- 
energy tensor of the sources supporting the Universe, one 
can prove that the Universe must be geodesically incom- 
plete ("singular"). Even in scenarios where the current 
ACDM cosmology was preceeded by a phase of slow- 
roll inflation [5], with eternal inflation occurring on even 
larger scales, it is a striking result t 6j that the initial sin- 
gularity remains. 

It is instructive to discuss which energy conditions need 
to be assumed to prove existence of a cosmological sin- 
gularity for the FLRW cosmologies 
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For k = —1,0 the only condition that must be assumed 
is the null energy condition (NEC), i.e. eqn. where 
w M is a future-pointing null vector field. The NEC is 
reasonable and in agreement with the known macroscopic 
matter and energy sources in our Universe |19j . 

For k = +1, however, the strong energy condition 
(SEC) must be assumed [20]. We are essentially certain 
that this condition is violated by macroscopic sources 
in our world, as well as in many completely consistent 



theoretical toy models. The goal of this paper is to ex- 
plore the two questions above for k — +1 Universes with 
sources satisfying the NEC (but violating the SEC). We 
will find that one can make classical cosmologies that live 
eternally, undergoing an infinite sequence of non-singular 
bounces, and remaining within the regime of validity of 
general relativity. When the ratio between maximal and 
minimal scale factors is not too large, these cosmologies 
are stable to small perturbations. In the opposite regime, 
when the ratio is large, we instead find both classical and 
quantum pathologies; classically there are growing modes 
(which can be tuned away), and quantum mechanically, 
particle production backreacts significantly after some 
number of cycles, likely causing a singular crunch 21j. 

Solutions. The FRW equations for the metric eqn. ^ 
are 
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where p is the energy density and p is the pressure. We 
want oscillatory solutions, namely those with two ex- 
trema (a = 0) such that at the smaller (where a = a_) 
a > 0, and at the larger (where a = a + ) a < 0. It is 
easy to see that these requirements, along with the NEC, 
only allow solutions for a when there is positive curva- 
ture, k = +1. The minimal model which oscillates has 
three components: positive curvature, a negative cosmo- 
logical constant (energy density = A < 0), and a "mat- 
ter" source with equation of state in the range 



p = wp, — 1 < w < — 1/3 



(4) 



(we will see a bit later that it is important that this 
source not be a perfect fluid). For this content the en- 
ergy density is p = A + p a a~ 3 ( 1+w ) where po is a constant 
parametrizing the density of the "matter." Then the so- 
lution to eqns. ^ is oscillatory. 

In the special case that w = — g these equations just 
describe a constrained simple harmonic oscillator and the 
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solution (setting k = +1) is 

a = + ao cos (wi + tp) 
where ip is an arbitrary phase and 



CO 



87r ^,A, 1 / 3A 

y G|A|, ao = _^_ 



(5) 



(6) 



This requires Po — ^ or P°sitivity of the radicand. 

Note that the Universe is static when this condition is 
saturated, though this requires a fine-tuning. In the op- 
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posite limit, Q — > oo, the ratio of the maximum to the 
minimum sizes a+ /a_ of the Universe goes to infinity. 

It is useful to switch to conformal time 77, where dn 2 = 
dt 2 /a(t) 2 . Defining 
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the solution for the scale factor ^ becomes 



a(rj) = - 



uj 1 — -^/l — 7Cos(?7) 
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Here lo is the frequency of oscillations given in ([5]), and 
we have set ip = 0. Notice that 7 ~ 4a_ /a+ for small 7. 

Stability. There are several simple stability issues we 
discuss here. (See e.g. [7] for a discussion of the corre- 
sponding stability issues in the Einstein static Universe.) 
First of all, the "matter" source in Eqn. Q may itself 
present dangers. In fact the canonical source which be- 
haves this way, a perfect fluid, would present a serious 
problem. To see this, recall that for scalar perturbations, 
one considers a more general metric 



ds 2 = a{ri) 2 [-(1 + 2$(r?, x))drj 2 + (1 - 2^(r], x))dQ%] 



For perfect fluids, $ = ^, Sp — c 2 5p, and 
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The derivatives are with respect to conformal time, and 
H = a! I a. An important point, clear from the sign of 



the V| 3 term in (10 1, is that if c 2 is negative, there are 
disastrous short-distance (high-momentum) instabilities. 

Now, a perfect fluid with w < —1/3 would have nega- 
tive c 2 . However, as explained in [8], one can find matter 
sources supporting equations of state of the form Q but 
with c 2 > (and in fact comparable to the speed of 
light), if one considers a "solid" with elastic resistance to 
shear deformations. A canonical example which they dis- 
cuss is a frustrated network of domain walls, which in the 
leading approximations gives precisely the w = — | case 
with the simplest behavior of a(t). For our purposes, 
the crucial point is simply that once we have achieved 



c 2 sufficiently positive, it is easy to check that the scalar 
perturbations above are stable. Further discussion of this 
point will appear in [9]. 

In addition to the above scalar perturbations, we need 
to consider tensor perturbations. These are governed 



by an equation whose form is identical to that of ( 13 ) 
below, and will be analyzed there. Next, homogenous 
but anisotropic perturbations are given by the Bianchi 
type IX metric [TD] ds 2 = -dt 2 + £- = i a 2 {t)a 2 , where 
are the Maurer-Cartan forms on S 3 . It is useful to 
parametrize the a, by an overall a(t) and two 'shape' 
deformations /3±(t), 
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Linearizing the FRW equations for f3± <C 1 then obtains 



01 + 2H(3' ± + 8k/3 ± = . 



(12) 



These modes will be analyzed momentarily. 

Another potential source present in our Universe, and 
indeed in any cosmological scenario, is gravity itself. 
That is, the Universe may respond to a produced gas 
of gravitons. The dynamics of massless particles may be 
described by a probe scalar field, with equation of motion 
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(13) 



We note as a curiosity that because of the periodicity 



of a, ( 10 1 and ( 13 1 can be recast as a Schrodinger prob- 



lem characterising motion of electrons in a particular Id 
periodic potential (where Bloch's theorem applies). 



The three types of perturbations (10 1, (12 1 and (13) 



have a similar structure; in fact, the anisotropic pertur- 



bation ( 12 1 is just a particular case of ( 13 ). Tensor modes 
of the metric are also described by eqn. (13). We denote 



a generic linearized mode by u, and expand in spherical 
harmonics, V^ 3 it; = — 1(1 + 2)u;. We now summarize the 
results of our numerical analysis of perturbations. 

There are three regimes of momenta where we expect 
(and shall find) different behavior. It is important to 
distinguish Universes with 7 ~ 0(\) from those with 
7 <C 1; we shall describe the behavior in both limits. 
• I = homogeneous mode: we expect that shifting such 
a mode should be analogous to shifting the homogeneous 
mode of the scale factor, which would simply move us 
in the space of periodic solutions and lead to a linear 
growth of the perturbation in naive perturbation theory 
(since e.g. two sinusoidal functions with slightly different 
frequency will perturbatively grow apart at a linear rate, 
as they get out of phase). This is borne out by the nu- 
merics for both 7 < 1 and 7 ~ 1. Thus, although this 
looks like a growing perturbation, that is likely just a 
failure of perturbation theory. Certainly a homogeneous, 
isotropic perturbation to the metric (as opposed to 4>) 
just moves us to a different one of our solutions and is 
not a dangerous instability. 
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• modes with momentum 2 < I < -4= on the 5 3 : these 
have long enough wavelength to detect the difference be- 
tween our cosmology and Minkowski space. For 7 ~ 1, 
i.e. a Universe which is "quivering" around a mean size, 
we find that they have oscillatory behavior and are thus 
stable. In contrast, for 7 <C 1, they can be unstable; we 
shall discuss bounds derived from their behavior below. 



• modes with I 3> -A=: these have small enough wave- 
length that they should barely detect the departures of 
our metric from flat space. As expected, they behave 
more or less like typical Minkowski space scalar field 
modes for times smaller than the period of oscillation 
of the Universe, for both 7 <C 1 and 7 ~ 1. 

The I = 1 mode is special. The perturbations gov- 
erned by (131 are stable for 7 ~ 1; on the other hand, 



the gravitational instabilities sourccd by ( 10 ) are always 
unstable for 1 = 1. For the case of a single-component 
perfect fluid, on which we have focused so far, this mode 
is absent from the physical spectrum: di^i = \ generates a 
global rotation on the S 3 and hence is pure gauge. How- 
ever, in multi-component systems there will generically 
be entropy perturbations; these contribute an inhomoge- 



neous term to ( 10 1 and can source a physical I = 1 mode. 



We find that the corresponding metric scalar mode ^1—1 
grows for all 7, unlike the case of modes with I > 2. [2"2"| 
However, we point out that even in these cases the 
I = 1 growing mode may be absent due to different 
mechanisms. A simple variant of our setup would be 
to orbifold the S 3 by a freely acting group in order to 
project out this mode. Orbifolding does not change the 
equations of motion (only local quantities appear there) 
but will project out modes from the spectrum. We have 
also not included the effects of non-gravitational damping 
modes. For example, a gas of gravitons will be subject to 
collisionless damping from free streaming, which damps 
growth at a rate proportional to the frequency Wfe of the 
mode. There is a range of 7 for which the growth rate of 



the I = 1 mode predicted by ( 10 1 is smaller than the cor- 



responding free streaming damping rate, thus killing this 
mode. The other fluids in the setup, including the do- 
main wall network, may also have other collisional forms 
of damping that can reduce the growth of this mode. We 
will discuss these points further in 0. In what follows 
we will assume that the / = 1 growing mode is absent. 

To summarize, the Universes with 7 ~ 1 are classically 
stable at the linearized level and live forever. The Uni- 
verses with 7 < 1 suffer from exponential growth (as a 
function of cycle number) of the finite momentum modes 
with I <C We show the numerical analysis of the 



modes of eqn. ( 13 ) in Figure[T]for all three regimes of mo- 
menta and various values of 7. The exponential growth 
whose beginning is shown in the middle figure would not 
be present for 7 ~ 1. The metric scalar perturbations 
\& behave in a qualitatively similar way, although they 
exhibit a faster rate of growth due to the gravitational 




FIG. 1: Massless scalar field evolution in conformal time, 
for different values of momenta. The first plot shows the 
homogeneous (I = 0) sofution with 7 = 10 . The second 
piot corresponds to I = 2 and 7 = 0.225; three eyefes are 
included, showing the exponential growth in the amplitude. 
The third plot has I = 45 and 7 = 0.01, and shows a single 
cycle. The initial conditions are (f)(0) = and (j>'(0) = 1. 



backreaction included in eqn. ( 10 ) 



Classical and quantum destruction of the Universe. 
For 7 ~ 1, the Universes we are studying are classically 
stable. For 7 <C 1, the exponential growth of the modes 
with < I < clearly indicates that we should expect 

such a Universe to have a bounded lifetime (at least un- 
til our approximations break down) . Can we tune this to 
allow a large number of oscillations within our regime of 
computational control? 

The cross-over from exponential to oscillatory behav- 
ior in the numerical solutions at I ~ l c = together 
with basic attempts to fit the growing solutions, suggest 
a rough form for the growing modes 



Ui(N) ~ u exp 




x N 



(14) 



where c ~ O(l), and ui(N) denotes the value of the Ith 
momentum mode after N oscillations, with starting vev 
uq. The important physical question is: when does the 
energy density in these modes become large enough that 
they compete with the dominant energy sources present 
in our background geometry? The ratio of the energy 
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density in the scalar perturbation to the cosmological 
constant is given by 
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Using ( 14 ) , and evaluating the resulting integral in a 



saddle-point approximation, we find the dominant I is 
Caddie ~ 'c/^> an d * ne ener gy ratio is thus 
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So, backreaction from the classical scalar field becomes 
important after a number of cycles N c given by 



log 



M P 
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(17) 



Classically, by tuning the initial state to have sufficiently 
small uq in Planck units, we can obtain an arbitrarily 
large lifetime even for the systems with 7 <C 1. 

Quantum mechanics is expected to induce an RMS 
value of uq, preventing a classical tune from saving the 
Universe for 7 < 1. Let us show this explicitly for the 



scalar field ( 13 ). To quantize the field, we impose canoni- 



cal commutation relations on the canonically normalized 
scalar x = a ( r t) c t > i 



bc(O),d v x(0')]= is {3 \e-e'), 



(18) 



where coordinatizes the three-sphere. This implies 
that in the instantaneous ground state characterizing the 
scalar at a time when the Universe has scale factor a, 



We are 



free to choose, as our initial quantum state, the instanta- 
neous vacuum associated to any value of the scale factor. 
Choosing, for instance, the "natural" quantum vacuum 
associated with a = a + (where the Universe is large and 
smooth and we have a natural expectation for the vac- 
uum state), gives 4>q ~ w^. This gives a bound on the 
number of cycles 



N c ~ log 



M P 



(19) 



This can be made parametrically large for small values 
of A (or, in this particular case, for sufficiently small 7). 
For 7 ru 1 the solutions to ( 13 ) are oscillatory, so the 



RMS values for various fields induced by quantum me- 
chanics will not cause instabilities. Hence for these values 
of 7, the universe is stable against perturbative classical 
and quantum instabilities. 

Conclusions and Questions. Our model with 7 ~ 1 
seems to provide an example of an eternal universe with- 
out singularities. It avoids the singularity theorems 
by having positive curvature and violating the SEC, 
though not the NEC. This universe is both classically 



and quantum mechanically stable against small pertur- 
bations at linearized level. Possibly, however, the back- 
ground "solid" could have microscopic dynamics that 
produce entropy, leading to a singularity even in our 
seemingly eternal models. This is an interesting, but 
model-dependent, question. We have focused on model- 
independent bounds here. This raises the question, can 
we find assumptions (weaker than the SEC) allowing 
proof of a 'quantum singularity theorem' that applies 
to closed Universes, extending the results of [3] to this 
physically important case without assuming unphysical 
energy conditions? We intend to pursue these and other 
questions in the future [9]. 



The cyclic nature of these cosmologies strongly sug- 
gests searching for exactly periodic quantum states in 
our geometry, perhaps characterizing special choices of 
the wavefunction of the Universe. Could some of these 
special quantum states be eternal, and provide "natu- 
ral" boundary conditions for certain closed cosmologies, 
in analogy with [llj? Perhaps this could even allow com- 
pactification of the time dimension. 



How realistic are these models? In particular, can 
we embed realistic ACDM cosmologies, with a preced- 
ing phase of inflation, into the expansion phase of one 
of our cycles in the 7 < 1 case? This would require a 
transition from radiation (and then matter) dominance 
during expansion to curvature (and our "solid") domi- 
nance near the following bounce. Given their relative 
scalings with a, this may require the radiation and mat- 
ter modes to be "Higgsed" above a large energy scale. 
As we have seen, such a universe with 7 <C 1 appears 
unstable. However, we were maximally pessimistic in ig- 
noring free streaming; could this effect vitiate the growth 
of inhomogeneous perturbations? Alternatively, for the 
stable, eternal 7 ~ 1 cosmologies, can we envision a Uni- 
verse which begins in such a phase, persists there for a 
long period, and then tunnels to a realistic inflationary 
Universe? Could either of these possibilities demonstrate 
that our observed universe might not have emerged from 
an initial singularity? 
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